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We report on the resummation of soft-gluon emissions for squark-antisquark production at next- 
to-next-to-leading-logarithmic (NNLL) accuracy. We will put particular emphasis on the one loop 
hard matching coefficients required to perform the resummation. Furthermore we will discuss the 
numerical effect of the different ingredients in the corrections. We find a significant reduction in 
the scale uncertainty and a considerable increase in the prediction of the total cross section at the 
central scale. Compared to the next-to-leading order prediction, the corrections increase the cross 
section by up to 30% for 1.5 TeV squarks at a centre-of-mass (CM) energy of 7 TeV at the LHC. 
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1. Introduction 

Supersymmetry (SUSY), and in particular the minimal supersymmetric Standard Model [1], 
can provide a solution to the hierarchy problem, accommodate gauge coupling unification and offer 
a dark matter candidate. Since these issues are addressed only if the SUSY scale is comparable to 
the weak scale, one would expect new SUSY particles (sparticles) with masses in the TeV range, 
which could be measured at the Large Hadron Collider (LHC). Particularly the coloured sparticles, 
squarks {q) and gluinos {§), would be pair-produced copiously in hadronic collisions and thus offer 
the strongest sensitivity. Searches at the LHC with a centre-of-mass (CM) energy of \/S = 7 TeV 
have placed lower limits on squark and gluino masses around 1 TeV [2, 3]. 

The leading order (LO) predictions for inclusive squark and gluino hadroproduction depend 
strongly on the renormalization and factorization scale. This dependence is reduced significantly if 
higher-order SUSY-QCD corrections are included. If the renormalization and factorization scales 
are chosen close to the average mass of the produced particles, the corrections generally increase the 
size of the cross section compared to the LO prediction. Consequently, the SUSY-QCD corrections 
have a substantial impact on the determination of mass exclusion limits and lead to a significant 
reduction of the uncertainties on SUSY mass or parameter values in the case of discovery [4,5]. 

The squark-antisquark production processes have been known for quite some time at next-to- 
leading order (NLO) in SUSY-QCD [6,7]. A significant part of the NLO QCD corrections can 
be attributed to the threshold region, where the partonic CM energy is close to the kinematic pro- 
duction threshold and the NLO corrections are dominated by soft-gluon emission off the coloured 
particles in the initial and final state. These soft-gluon corrections can be taken into account to all 
orders in perturbation theory by means of threshold resummation techniques [8,9]. This has been 
done for all MSSM squark and gluino production processes at next-to-leading-logarithmic (NLL) 
accuracy [10-14]. 

Recently, resummation at NNLL accuracy for squark-antisquark pair production was pre- 
sented [15]. Compared to the NLL calculation the new ingredients are the one-loop matching 
coefficients, which contain the NLO cross section near threshold, and the two-loop soft anomalous 
dimensions. We wiU discuss the impact of the corrections and provide an estimate of the theoretical 
uncertainty due to scale variation. 

In section ^ we discuss NNLL resummation for squark-antisquark pair-production. In sec- 
tion H we present the calculation of the hard matching coefficients. The numerical results for the 
LHC with a CM energy of = 7 TeV are presented in section ^ We conclude in section ^. 

2. Threshold resummation at NNLL 

We first briefly review the formalism of threshold resummation for qq production. The inclu- 
sive cross section C7/, can be written in terms of its partonic version C/y^^i as: 

(^hh2^q-q{pAfn^}) = L / dxidx2dpd(p - — ) fi/h,ixi,n^)fj/h2ix2,l^^)(yij^^^{p,{m^},n^) , 

j j J V XiX2 J 

where {m^} denotes all masses entering the calculations, /,7 are the initial parton flavours, /,y/j| 
and /y//,, the parton distribution functions, \i is the common factorization and renormalization scale 
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and x\ 2 are the momentum fractions of the partons in hadrons h\ and h2. The hadronic threshold 
for the inclusive production of two squarks with mass corresponds to a hadronic CM energy 
squared S = 4m^. Therefore we define a threshold variable p = 4m?/5 and its partonic equivalent 
p = Aniq/s, with s = x^xiS the partonic CM energy squared. The resummation of the soft-gluon 
contributions is performed after taking a Mellin transform (indicated by a tilde) of the cross section, 



(2.1) 



The threshold limit p — )• 1 corresponds to A/^ — )• oo. Near threshold, fixed-order perturbation theory 
does not converge well due to logarithmically enhanced terms. By resumming large logarithmic 
corrections containing L = \og{N) to all orders, a new perturbative expansion arises, where first the 
leading logarithmic (LL) corrections are resummed, followed by the NLL, NNLL, ... contributions. 
The relation between this new perturbative expansion and fixed-order perturbation theory is shown 
schematically in Figure |]. 
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Figure 1: Schematic form of the logarithmically enhanced terms with the LO cross section factored out. 
The terms enclosed in the solid line occur in the LL approximation, which completes the first column. The 
new terms in the NLL approximation are enclosed in the dotted line and complete the second column. The 
shaded region corresponds to new terms coming from the NNLL approximation, which completes the third 
and fourth column. We resum the logarithms to all orders, so the columns extend downwards to infinity. 

The all-order summation of these logarithmic terms is based on the near-threshold factorization 
of different classes of radiation: hard, (soft)-collinear, and wide-angle soft radiation [8,9, 16-19]. 
Near threshold the resummed partonic cross section to NNLL accuracy has the form [8,9,20]: 



JO) 



ij^qqj 



1 + 



«s ^^Coul,(l) |'.„2t ,,2- 



j{N,{m^},H^) exp Lgiia^L) + gjjia^L) + asg3,iia,L) 



K 



1 + ^.^(1) ({^2| 2 



(2.2) 



where O'.^j^ -= ^ are the colour-decomposed LO cross sections in Mellin-moment space. The colour 
label / corresponds to an irreducible representation of the colour structure of the process, which for 
squark-antisquark production can be either a singlet or an octet [10, 11,21]. The exponent in the 
first line of Eq. ( |2.2[ ) captures all dependence on the large logarithm L. The last line contains the 
one-loop Coulomb contribution ^*^°"''(^) and hard matching coefficient 

Setting the hard matching and Coulomb coefficients to and keeping only the g\ term in 



Eq. (2.2) constitutes the LL approximation. Including the g2 term as well corresponds to NLL. 
For NNLL accuracy also the g3 term and the one-loop hard matching and Coulomb coefficients 



need to be taken into account, as can be seen explicitly by expanding Eq. (2.2) and comparing it to 
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figure H. The §3 term can be found in [15, 22-24] and the ingredients for in [15, 25, 26]. 

The calculation of the NLO hard matching coefficients ^.^.^^ -= , will be discussed in section |3l 

The hadronic cross section Of,^ij^^^^j can be recovered by taking the inverse Mellin transform. 
We retain the information contained in the complete NLO cross section [7] by combining the NLO 
and NNLL results through a matching procedure that avoids double counting of the NLO terms: 

<;°;-r""''""'('^'{-'}'^') =1 L 'iiP-'f'/'niN+l,^')fjM^+l,^^ (2.3) 

We adopt the "minimal prescription" of Ref. [27] for the contour CT of the inverse Melhn transform 
in Eq. In order to use standard parametrizations of parton distribution functions in ;c-space 

we employ the method introduced in Ref. [28]. 

3. Calculation of the hard matching coefficients 

We now discuss the calculation of the one-loop hard matching coefficients The calcula- 
tions were done using FORM [29]. After performing an expansion of the NLO cross section in jS, 
the hard matching coefficients ^'^^ are determined by the terms in the NLO cross section that are 
proportional to j3, j8log(j8) and /3log^(j8). Terms that contain higher powers of j5 are suppressed 
by powers of l/N in Mellin-moment space and do not contribute to the matching coefficient. In 
contrast to the case of top-pair production in Ref. [30], there is no full analytic result for the real 
corrections to squark-antisquark production, so we cannot take the explicit threshold limit. 

The virtual corrections for squark-antisquark production can be obtained from the full analytic 
calculation as presented in Ref. [7]. First we need to colour-decompose the result. We only need 
the colour decomposition of the LO matrix element. Due to the orthogonality of the ^-channel 
colour basis that we use, the full matrix element squared is then automatically colour-decomposed: 



NLO,/ — 2Re(^NLO-^LO,/ 



Then we are left with an expression in terms of masses, Mandelstam variables and scalar integrals. 
Since we need the cross section to ^(j8), we have to expand to zeroth order in jS to obtain 
the virtual part of the hard matching coefficients. 

The integrated real corrections at threshold are formally phase-space suppressed near threshold 
unless the integrand compensates this suppression. Therefore we can construct the real corrections 
at threshold from the singular behaviour of the matrix element squared, which can be obtained 
using dipole subtraction [31,32]. First we recall that the cross section can be split into three parts: 
a part with three-particle kinematics o^^\ one with two-particle kinematics o^^\ and a collinear 
counterterm a'' to remove the initial-state collinear singularities. These parts are well-defined in 
n =4 — 2e dimensions, but their constituents diverge for £ — 0. With the aid of an auxiliary cross 
section a^, which captures all singular behaviour, all parts are made finite and integrable in four 
space-time dimensions. This auxiliary cross section is subtracted from the real corrections at 
the integrand level to obtain o^^^ and added to the virtual corrections a^, which defines o^^^: 
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Compared to the case of two-parton kinematics, the phase space of a^^^ is limited by the energy 
of the third, radiated massless particle. Near the two-particle threshold, the maximum energy of 
the radiated particle, and thus the available phase space, equals ^max = \fs — lmq oc j8^. Since after 
subtracting no divergences are left in the integrand of a^^J^, the leading contribution of o^^^ is at 
most proportional to and can thus be neglected. So at threshold the real radiation is completely 
specified by the singular behaviour contained in and can therefore be determined by summing 
over dipoles that correspond to pairs of ordered partons [32] and taking the threshold limit. 

After combining the real and the virtual corrections the hard matching coefficients can be ob- 
tained by taking the Mellin transform and omitting the Coulomb corrections and the log [N) terms. 
The complete expressions for the hard matching coefficients of the squark-antisquark production 
processes can be found in Ref. [15]. Their behaviour for varying gluino mass is shown in Fig. ^. 



0.7 I : ^ : ^ 1 0.5 




500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 

(a) "'9 (b) '"9 

Figure 2: Gluino-mass dependence of the colour-decomposed NLO hard matching coefficients for the qq 
initiated channel (a) and the gg initiated channel (b). We have set fi = nig = 1.2 TeV and nit ~ 172.9 GeV. 



4. Numerical results 

We now present numerical results for NNLL-resummed squark-antisquark pair-production at 
the LHC for a CM energy of 7 TeV. We use the 2008 NLO MSTW parton distribution functions [33] 
with the corresponding as{M^) = 0.120. We have used a top quark mass of nif = 172.9 GeV [34]. 
In order to study the effects from the hard matching coefficients and the Coulomb corrections sep- 
arately, we compare the NLL matched cross section (jNLO+nll^ which is based on the calculations 
presented in [10-12], the NNLL matched cross section without the Coulomb contributions to the 
resummation a^LO+NNLLw/o Coulomb^ ^Yiich has <rc°"''(i) = in Eq. Q, and the NNLL matched 
cross section (jNLO+nnll^ which does include the Coulomb correction <^'^°"''(i). All cross sections 



are matched to the NLO resuh calculated with PROSP INO [7, 35] using Eq. 

Figure |3|(a) shows the mass dependence of the scale uncertainty for the squark-antisquark cross 
section. The squark and gluino mass have been taken equal and the scale has been varied in the 
range mq/2< pL < 2m q. We see that the scale uncertainty reduces when higher-order corrections are 
included. In this range of squark masses, the NNLL resummation without the Coulomb corrections 
<^Coui,(i) already reduces the scale uncertainty to at most 10%. The inclusion of the Coulomb term 
<^Coui.(i) (.jjg resummed NNLL prediction results in a scale uncertainty of only a few percent. 
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Figure 3: The scale uncertainty (a) and K-factor (b) of the LO, NLO, NLO+NLL and NLO+NNLL (both 
with the Coulomb part '^'^o'^' ti) and without it) squark-antisquark cross sections for the LHC at 7 TeV. The 
central scale and the gluino mass have been taken equal to the squark mass. In (a) the common renormaliza- 
tion and factorization scale has been varied in the range < |i < 2mq. 



In figure ^b) we study the mass dependence of the .^-factor Kx = o^/ a^^*^ at the central scale 
H =mq and rrig = mq. The theoretical prediction of the cross section increases as more corrections 
are included. The effect is more pronounced for higher masses, since in that case particles are 
produced closer to threshold. The NNLL resummation without the Coulomb corrections •^'^""' (i) 
already results in a 25% increase of the cross section with respect to the NLO cross section for 
squarks of 2 TeV. The contribution from the Coulomb term to the resummed NLO+NNLL cross 
section is even larger, yielding a total K-ia.ctor of 1.45. Although the effect from the Coulomb 
corrections could be somewhat smaller in reality due to the finite lifetime of the squarks, figure ^b) 
suggests that the NNLL contribution will remain large. 



5. Conclusions 

We have discussed the NNLL resummation of threshold corrections for squark-antisquark 
hadroproduction and presented numerical results for the NLO-i-NNLL resummed cross section for 
squark-antisquark production at the LHC with a CM energy of 7 TeV. For a squark mass of 2 TeV, 
the NLO-i-NNLL squark-antisquark cross section is up to 45% larger than the corresponding NLO 
cross section. The correction is reduced to 25% if the resummation contributions due to Coulombic 
interactions are omitted. The scale dependence is reduced significantly, particularly after inclusion 
of the Coulomb corrections. However, the observed reduction in the scale dependence might be 
modified somewhat by the inclusion of the width of the particles or by matching to the full NNLO 
result, which is not available. A very conservative estimate of the scale uncertainty is provided 
by the NLO-i-NNLL w/o Coulomb results, which do not include the Coulomb corrections. This 
information could be used to improve current limits on SUSY masses or, in the case that SUSY is 
found, to more accurately determine the masses of the sparticles. 
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